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ABSTRACT 

We show that invaluable information on the structure quasar outflows can be obtained by considering mi- 
crolensing induced variability of absorption line troughs in lensed quasars. Depending on the structure and 
geometry of the outflowing gas, such extrinsic line variability can be manifested as changes to the equivalent 
width of the line as well as line profile distortions. Here we consider several physically distinct outflow models, 
having very similar spectral predictions, and show how ML induced absorption line variability can be used to 
distinguish between them. Merits of future systematic studies of these effects are exemplified. 

Subject headings: gravitational lensing — ISM: jets and outflows — galaxies: active — galaxies: nuclei — 
quasars: absorption lines 



1. INTRODUCTION 

A large fraction of type-I (broad emission line) quasars 
show gas outflowing from their centers with velocities in 
the range 10 2 — 10 kms~' (inferred from the detection 
of blueshifted absorption lines; see Crenshaw, Kraemer, & 
George 2004 and references therein). Such flows are believed 
to be launched from the central engine of quasars (e.g., Elvis 
2000) and may have a considerable effect on their environ- 
ment (e.g., King 2003, Scannapieco & Ho 2004). 

Physical models for quasar outflows include radiation pres- 
sure driven cloud/wind models (e.g., Arav, Li, & Begelman 
1994, Murray, Chiang, Grossman, & Voit 1995, Proga, Stone, 
& Kallman 2000, and Chelouche & Netzer 2001), magnetic 
winds (e.g., Konigl, & Kartje 1994), and thermal pressure 
driven flows (e.g., Woods et al. 1996, Chelouche & Netzer 
2005). While different models imply very different physical 
properties of the outflowing gas (such as location, density, and 
mass loss rate), they are all in qualitative agreement with ob- 
servations. Thus, despite some 30 years of research little is 
known with confidence about the physics of such flows. 

In Chelouche (2003) we have shown that the effect of strong 
gravitational lensing can be used as a probe for the structure 
of quasar outflows over angular scales comparable to the sep- 
aration between quasar images. The proposed method was 
useful in cases where the effects of gravitational microlensing 
(ML) may be neglected. The focus of this work is to consider 
cases where ML effects are likely to dominate (e.g., Ofek & 
Maoz 2003) and investigate how they can be used to shed light 
on the physics of quasar outflows. Specifically, we focus on 
the information that can be gained by studying the effects of 
ML on absorption line profiles in lensed quasars. Such effects 
are potentially significant (e.g., Hutsemekers 1993, Lewis & 
Belle 1998) yet the implications for studying quasar outflows 
in the framework of realistic physical models have not been 
addressed. 

This paper is organized as follows: In §2 we discuss current 
paradigms in flow modeling and elaborate on the formation 
of absorption line profiles in quasars. Section 3 describes the 
effect of ML on the shape of line troughs for several flow 
models. Conclusions follow in §4. 
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2. THE FORMATION OF LINE TROUGHS 

The formation of absorption line profiles depends on the 
projected 2D properties of the absorbing gas (such as velocity, 
temperature, and ionization structure) and on the continuum 
source (e.g., emissivity and size). The residual flux under the 
line formed by a flow which is located between us and the 
continuum emitting source satisfies, 



■Fobs(v) 



■■ J dR 2 n(R)I{R)e 



-t(R.i 



(1) 



where R is a vector marking the projected position over the 
plane of the sky, / the emitted flux per unit projected area and 
is a slowly varying function of wavelength in the vicinity of 
the line (for the sake of simplicity we neglect the possible con- 
tribution from emission lines), jti(R) the (de-)magnification 
pattern due to ML by stars in the lensing galaxy (see below), 
and t(R,v) the position and velocity dependent optical depth 
which depends on the properties of the outflowing gas. 

An equivalent, observationally inclined expression for the 
residual flux is 



F obs (v)=F{[l-C(v)]+C(v)e-^}, 
where the unabsorbed flux emitted by the source, 
F= f dR 2 I(R). 



(2) 



(3) 



t(v) and C(v) are the observed velocity dependent optical 
depth and covering factor, respectively (e.g., Arav, Korista, 
& de Kool 2002). A major breakthrough in flow research was 
the realization that the shape of the rest ultraviolet (UV) line 
troughs (e.g., C IV 1549A) is determined by velocity depen- 
dent partial covering effects (e.g., Arav et al. 1999; hereafter 
A99). Both t(v) and C(v) can, in principal, be deduced from 
observations and used to infer on the physics of the flows (i.e., 
on t(R, v), A99). In practice, this inversion problem is under- 
constrained and physically distinct outflow models give simi- 
lar line profiles. This model degeneracy is explained below. 

2.1. Ordered flows 

This model, as proposed by A99, assumes that the flow 
is inclined to our line of sight by angle 9 and is monotoni- 
cally accelerated across the projected surface of the contin- 
uum source (see the right panel of figure 1). This means that 
every projected velocity interval, 8v (= 5vcos0, where v is 
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accretion disk 

FIG. 1. — Ordered (A99) vs. disordered (de Kool et al. 2002) flow geometries. Ordered flows (right panel) are envisioned as monotonically accelerating 
gas streams (with thickness /) that cross our line-of-sight with inclination angle 9. Disordered flows (left panel) may consist of numerous clouds randomly 
distributed in real space and in velocity space, each having its properties (e.g., optical depth) drawn from the same parent population. We study cases where the 
total line-of-sight covering factor of the flow, Cn ow , is equal or smaller than unity (see §3.1). 
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FIG. 2. — The absorption line profiles predicted by ordered and disor- 
dered flows. Note the similarity between the lie troughs predicted by the 
different models. Here we assume: Cfl ow = 1, v c = 1200 km s _1 , C max = 
0.63, Av = 5000 km s _1 , T 3> 1 for the ordered flows and p = -0.62 and 
T m „(v< 1200kms"Vv 2J ; T max (v > 1200 km s _1 & v < 3800 km s _1 ) = 
10 T m ax(v > 3800 km s _1 ) « y~ 2 - 3 for the disordered flows. 



measured along flow lines; see figure 1), corresponds to a dif- 
ferent region in the source plane. (Note, however, that in this 
model there is a range of velocities at every location R). Here 
we consider a simplified version of the model where the ve- 
locity gradient, dv/dr (r is measured along flow lines) and the 
flow density are constant. For a flow whose line troughs ex- 
tend from v = km s _1 (corresponding to the rest wavelength 
of the line) to some maximum velocity, Av, we obtain 

t(v) = const. 

{C m axV/v c ; V<V C 
C max ; v > v c and Av - v > v c 

Cmax (Av - v)/v c ; Av - V < v c 

(4) 



where C max = IcosO/R^ is the maximum covering factor, and 
v c = (dv/dr)(l /tg0)cos0 . Cflow is the fraction of the ionizing 
source which is covered by the outflow. An example of the 
absorption line profile produced by such a flow is shown in 
figure 2. 

2.2. Disordered flows 

A different model, suggested by de Kool et al. (2002; cf. 
Lewis & Belle 1998), considers a flow comprised of numer- 
ous small (compared to Rd) clouds that are randomly dis- 
tributed in real and in velocity space. Each cloud from an 
ensemble of clouds traveling with velocity v (there are many 
such ensembles covering the full velocity range, Av) has a dif- 
ferent optical depth drawn from a powerlaw distribution with 
index p and a maximum optical depth distribution, T max (v) 
(see equation 9 in de Kool et al. 2002). The residual flux (see 
their equation 13) is 

F obs (v) = F ^ + l ™ + l] P [p+l, W (v)] • (5) 

T max (v), p are not known a-priori, but can be deduced from 
observations, and F, P are the complete and incomplete 
gamma functions, respectively (see de Kool et al. 2002 for 
a detailed layout of the formalism). Figure 2 demonstrates 
that, the absorption line profiles from ordered and disordered 
flows can be very similar for a suitable choice of parameters, 
and cannot therefore be easily distinguished. Below we wish 
to examine whether this degeneracy can be removed by con- 
sidering the effects ML has on such line profiles in lensed 
quasars (e.g., Hutsemekers 1993). 

3. DISTORTION OF LINE TROUGHS BY MICROLENSING 

Recent studies have shown that ML can induce significant 
(> 0.1 magnitude) fluctuations in the light-curves of quasars 
(e.g., Irwin et al. 1989, Wozniak et al. 2000, Ofek & Maoz 
2003). This effect results from the Einstein radius of stars in 
the lensing galaxy being comparable to the size of the quasar 
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FIG. 3. — The magnification map, /i(R), used here to calculate the dis- 
tortion of the absorption line profiles formed by quasar flows (equation 1). 
Dark shades correspond to higher magnifications. A web of caustics is clearly 
shown resulting from the combined effect of many stars in the lensing galaxy. 
Here we assume a microlensing mass surface density, K = 0.55 with no exter- 
nal shear (y = 0). The rectangles mark the initial position and size, Rj, of the 
quasar continuum emitting region. Rj = 1 corresponds to roughly three times 
the Einstein radius of one solar mass star (~ 3 x 1() 16 cm for a typical lensed 
system where the quasar, lens redshifts are 2, 0.5, respectively; cf. Treyer & 
Wambsganss 2004). 



continuum source (e.g., Schneider, Ehlers, & Falco 1992 and 
references therein). It has been suggested that ML can be used 
to explore the central engines (e.g., Agol & Krolik 2000) and 
the broad emission line regions in quasars (e.g., Abajas et al. 
2002, Chartas et al. 2004, Lewis & Ibata 2004, Richards et 
al. 2004). In particular, Hutsemekers (1993) and Lewis & 
Belle (1998) have shown that single-caustic crossing events 
can distort the shape of absorption line profiles. Such dis- 
tortions result when the magnification pattern is non-uniform 
across the continuum source (i.e., ju(R) ^const in equation[0 
and can be manifested as equivalent width (EW) changes (cf. 
Lewis & Ibata 2003), absorption line profile distortions, and 
their combination. 

Qualitatively, pure EW changes occur when a continuum 
region which is occulted by the absorbing gas (e.g., a cloud) 
is (de-)magnified with respect to the non-occulted part of 
the source. Shape distortions occur when different parts of 
the flow traveling with distinct velocities cover separate re- 
gions of the continuum source which are subjected to dif- 
ferent ML magnifications. A more quantitative treatment re- 
quires the calculation of realistic magnification patterns and 
flow models (e.g., Kochanek 2004). The magnification map 
(source plane), ju(R) (see figure 3) was calculated using a 
ray-shooting algorithm similar to that which is described in 
Wambsganss (1999). Here consider a canonical case where 
the redshift of the lensed quasar, z s = 2, and the lensing galaxy 
is at, zi — 0.5 (such values are typical of strongly lensed sys- 
tems; see the CfA-Arizona Space Telescope Lens Survey 3 ; 
Munoz et al. 1998). We assume that the dimensionless sur- 
face mass density, K = 0.55 (see Schneider, Ehlers, & Falco 
1992), comprised of equal mass stars whose Einstein radius is 
~3x 10 16 cm for the chosen parameters (see figure 3). We 
assume no external shear. 

As we shall see below, it is instructive to consider three 
length scales of the problem: The size of the emission re- 

3 See http://cfa-harvard.edu/castles 



gion, R c j, the auto-correlation length scale of the magnifica- 
tion pattern, R^, and the auto-correlation length scale of the 
flow velocity field, R v . The latter scales are derived from the 
autocorrelation function, 

A x (r)oc J dr'X(r-r')X{r'), (6) 

where X is the property in question, and Rx is defined as the 
point beyond which Ax drops by an e-fold. Here r is measured 
in the mean (over the surface of the emission region) direction 
of the velocity gradient of the flow. Below we consider sim- 
plified models pertaining to the simplified versions of ordered 
and disordered flows and calculate the effect of ML on the 
absorption line profiles. We consider cases where ML effects 
are likely to be significant (e.g., producing detectable varia- 
tions in the light-curve of the quasar), i.e., where R^ ~ Rj. 

3.1. Ordered flows 

Here we consider cases which correspond to ordered flows 
(R v ~ R c /). We first assume that the flow completely covers 
the continuum source (Cflow = 1)- The effect of ML on the 
absorption line profiles formed by such flows is shown in fig- 
ure 4. Our primary conclusion is that line profile distortions 
can be considerable (> 10% flux changes; see also §3.4) and 
show a clear distortion pattern extending over the full velocity 
range covered by the flow. This pattern reflects the correlation 
between the flow velocity and the magnification pattern (both 
of which are correlated with R). This result is very different 
from that of Lewis & Belle (1998) and from that which cor- 
responds to disordered flows (see below) where line profile 
distortions occur over small velocity intervals. It is also noted 
that, for this flow configuration, no significant ML induced 
EW changes take place. This property is common to all lines 
that are formed by flows for which t(R, v) = t(v), i.e., whose 
properties (such as temperature or ionization) do not depend 
significantly on R (see appendix). 

The magnitude and shape of the distortions depend on the 
flow properties, t(v), C(v). Specifically, larger line optical 
depths correspond to larger distortion amplitudes (10% — 30% 
for t = 0.5 — 3 respectively, see figure 4a). This is also the 
case for larger covering fractions (shown in figure 4a as a 
function of C max )- This results from the fact that the lower 
residual flux in such cases is emitted by an overall smaller re- 
gion of the continuum source. Such regions are likely to expe- 
rience larger amplitude (de-)amplifications. Conversly, larger 
continuum emitting regions (or denser caustic maps due to a 
higher stellar density in the lensing galaxy; e.g., Wambsganss 
1992) result in lower amplitude, more complex line distor- 
tions (see figure 3b where different source sizes where used 
in each case). 

In terms of Rd,R^, and R v this behavior can be understood 
as follows: When R^ is of the order of R<] (~ R v , for ordered 
flows), large (de-)magnifications and considerable line distor- 
tions and EW variations are expected. For the specific exam- 
ple considered here, R^ ~ 0.5 (in units of Rd) 4 - (Note that R^ 
is sensitive to the value of K and may vary by a large factor 
for K in the range 0.2 — 0.8; see e.g., figure 1 in Treyer & 
Wambsganss 2004). For R^ <C Rj (see for example figure 4b 
where a source of size R ( / = 3 was assumed), different parts 
of the continuum emitting region are (de-)magnified by un- 
correlated magnification patterns so that line variations would 

4 Care must be taken when calculating Sjj since the clumping of caustics 
in the magnification pattern is a highly non-linear effect. Here R^ is averaged 
over the entire magnification map 
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have a more complex, multi-humped pattern and low resolu- 
tion spectroscopy is likely to find an overall suppression of the 
effect. On the other extreme where ^> R c /, line distortions 
would also be suppressed as different parts of the continuum 
source are magnified by more or less the same factor (though 
long-term ML induced light-curve fluctuations may still be 
considerable; see also §3.4). 

Until now we have considered cases where Cflow = 1. Nev- 
ertheless, flows may not fully cover the entire emission re- 
gion (e.g., Chelouche & Netzer 2005 and references therein. 
In such cases (Cflow < 1 and, effectively, t(R) ^const.) con- 
siderable EW changes may occur, as demonstrated in the top- 
left panel of figure 5 for the case where Cfl ow = 0.7. These 
changes would be statistically maximized if the non-occulted 
section of the continuum emitting region is continuous and 
~ R c [. We further discuss this issue in §3.3. 

3.2. Disordered flows 

Here we consider the distortion of absorption line profiles 
formed by disordered flows for which R v <C R c i- We first con- 
sider cases where Cfl ow = 1 and model those in the follow- 
ing way: An ordered flow is first constructed using the A99 
formalism. The flow is then divided into sections covering 
each a velocity interval 8v across (there are Av/8v ~ R,i/Rv 
such sections). A disordered flow is then obtained by ran- 
domly shifting the positions of its sections. The shape of the 
line profile is conserved by this shuffling procedure. Thus, R v 
is essentially the size of individual sections so that R v /Rd — 
8v/Av. For example, a model with R v /R c i ~ 5v/Av = 0.3 can 
be envisioned as a three stream flow similar to the one con- 
ceived by A99 (see their figure 10). 

Figure 4c demonstrates that the line profiles formed by less 
ordered flows have smaller distortion amplitudes and that the 
distortion pattern is more complex, having many low ampli- 
tude humps. This results from the fact that the velocity is 
poorly correlated with R and, hence, with the magnification 
pattern, For R v <C line distortions would be negligible 
even though the total ML (de-)magnification of the contin- 
uum source is large and continuum fluctuations substantial. 
In cases where R v ~ R^ <C Rd, the effect of gravitational mi- 
crolensing on the quasar lightcurve as well as on the line pro- 
files would be negligible. 

We note that although considerable line distortions are un- 
likely to occur for highly disordered flows, it is still possi- 
ble for such lines to exhibit EW variations if a considerable 
and continuous region of the emitting source remains non- 
occulted by the flow (Cfl ow < 1). EW changes would then 
be qualitatively similar to those shown in figure 5 although no 
ML induced line profile distortions would be observed in such 
cases. 

3.3. Time variability 

ML induced effects in the lightcurves of quasars may vary 
considerably with time (e.g., Schneider, Ehlers, & Falco 1992, 
Ofek & Maoz 2003). This results from the motions of the mi- 
crolensing stars within the potential well of the lensing galaxy, 
and the relative peculiar velocity between the quasar, the lens, 
and the observer. The relevant timescale for such variations 
is, roughly, the caustic-crossing timescale which depends on 
Rd and the peculiar velocity of the galaxy/stars. Typical 
timescales are of order a few months to years (e.g., Schnei- 
der, Ehlers, & Falco 1992) and are therefore longer then typ- 
ical dynamical timescales (e.g., Chelouche 2003). Here we 
wish to examine how time-dependent ML effects reflect on 



the observed absorption line profiles. We model such time 
evolution by artificially shifting the continuum emitting re- 
gion (with Rd = 1) across the magnification map (from left to 
right; see figure 3) and assume a stationary flow. 

Examples for the time-dependent variations of the line pro- 
files from ordered flows are shown for full (Cflow = 1) and 
partial coverage (Cfl ow = 0.7) models. Figure 5 demonstrates 
that line distortions evolve with time and a pattern is revealed 
whereby the dip at low velocities gradually shifts to higher ve- 
locities as the central caustic pattern (figure 3) sweeps across 
the source. As expected, the EW of the lines remains unaf- 
fected by line distortions for a full coverage of the source yet 
considerable EW changes do occur for the partial covering 
model (compare panels 1 and 5 in figure 4). 

Significant EW variations of the line troughs can also occur 
for the case of highly disordered flows provided a continuous 
section of the emission region remains non-occulted. We note 
however that in such a case, significant line distortions would 
be absent. 

3.4. Line profile variation probability 

ML induced distortions of absorption line profiles are of a 
statistical nature related to the position of the quasar over the 
magnification map. An interesting quantity is therefore the 
probability to observe some variation at any given time. To 
calculate this, we have randomly positioned the quasar over 
the magnification map (assuming a Rd = 1 emission region) 
and calculated the line variation for each case. We define the 
following quantities: The mean distortion amplitude over the 
line profile, the (local) maximum distortion amplitude, and 
the relative EW change (see the caption of figure [6] for the 
proper definitions). 

We first consider ordered flows where we find consider- 
able variations of the line profiles most of the time. For 
the case considered here, the mean distortion amplitude ex- 
ceeds ~ 10% most of the time (with mean distortions exceed- 
ing ~ 15% some 20% of the time). Maximum distortions 
of the line profiles over localized velocity intervals exceed 
20% most of the time and may be as large as 40% in some 
cases. Qualitatively similar results are obtained when assum- 
ing partial line-of-sight coverage although here considerable 
EW changes may occur (for the example shown here > 10% 
variations are frequent). 

Very low distortion amplitudes are expected for highly dis- 
ordered flows. Nevertheless, EW variations will occur if a 
continuous region region of the continuum source is not oc- 
culted by the flow (the result would be similar to that shown 
in figure 5 for the same covering factor). 

Our calculations show that there is no significant correlation 
between the continuum magnification (i.e., its light-curve) 
and the variation amplitude of the line profile. This is not 
unexpected given the fact that line variations require the mag- 
nification pattern to vary across the emission region while the 
total magnification does not (in fact, a coherent magnification 
or de-magnification pattern would produce the largest contin- 
uum variations). 

4. CONCLUSIONS & SUMMARY 

We have shown that the effect of ML on the absorption line 
profiles can be utilized to investigate the structure and dynam- 
ics of quasar outflows and remove the degeneracy between 
different outflow models when only one dimensional spectra 
is considered. Nevertheless, this requires the identification 
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FIG. 4. — ML induced line distortions shown as the ratio of the distorted profile to the un-distorted one (Cfl ow = 1 is assumed), (a) More optically thick lines 
(solid lines) or those with a larger covering factor (dashed lines) are more distorted by ML effects since the the residual flux is emitted by more localized regions 
of the source, (b) Larger sources (Rj 2> i?^) will be microlensed by a more complex and less correlated caustic pattern and the resulting line distortions are 
less prominent [sources of size = 1,2,3 (in fiducial units) correspond to magnitude fluctuations, Am = 0.3,0.1,0.03, respectively], (c) Line profiles of highly 
disordered flows (R v /Rd <C 1; see text) are less distorted since the magnification and velocity are poorly correlated. 
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FIG. 5. — Evolution of the line profiles for an ordered flow (assuming an Rj = 1 continuum source). Solid (dashed) lines assume full (70%) coverage of the 
continuum source (see text). Considerable line distortions are evident in both cases. Significant EW variations are expected for a flow that continuously covers 
only part of the continuum source (compare panels 1 and 5). In this example, a hump in the line profile seems to propagate with time from high to low velocities 
as a web of caustics (shown at the center of figure 3) sweeps across the source. 



of ML induced line distortion events which may be compli- 
cated by the fact that the line profiles of different quasar im- 
ages may differ due to reasons related to the flow (e.g., Che- 
louche 2003, Green & Aldcroft 2003, Ellison et al. 2004) or 
due to incorrect subtraction of light contribution from other 
sources (e.g., the lens). Spatially resolved spectral monitoring 
of lensed quasars is required to elevate some of the difficulties 
in identifying ML induced line variations. 

In cases where line variations are predominantly due to ML 
then the implications of the various distortion patterns for the 
structure and geometry of quasar outflows can be summarized 
as follows: 

• Significant line EW variations are observed: A local- 
ized, continuous region of the continuum source is not 
covered by the outflowing gas. 

• Significant line profile distortions are observed and 



those are well correlated with velocity: The flow cov- 
ers most of the continuum source and has a relatively 
ordered velocity structure. 

• Low amplitude, complex distortion pattern is observed: 
The flow is either disordered or the magnification pat- 
tern highly complex (in the latter case ML induced 
light-curve fluctuations would be suppressed too; e.g., 
Wambsganss 1992). 

Most promising candidates for observing line profile distor- 
tions are lensed quasars showing considerable extrinsic flux 
variations (i.e., R c j ~ R^) and having non-black, saturated line 
troughs. Of particular interest are those quasars having low 
velocity flows in which case t(v) and C(v) can be indepen- 
dently deduced (e.g., Arav, Korista, & de Kool 2002) and used 
to put further constraints on models. Spatially resolved spec- 
troscopy of two or more quasar images is required to distin- 
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0.2 0.3 
line profile variation 

FIG. 6. — The variation probability of the absorption line profiles for a 
given flow model and magnification map (a source of size one Rj is as- 
sumed). Three variation measures are defined: The mean variation is de- 
fined as (|0(v) [0(v) _1 - l] ^ where 0(v), 0(v) are the non-distorted, 
distorted line profiles, respectively. The maximum variation is defined as 
max j|0(v) [0(v)~' - l] | J. The EW variation is simply the relative EW 

difference between the ML and the non-ML troughs. Solid lines correspond 
to a model with Cfl ow = 1 (we note that no significant EW variations occur 
in this case). Large 10% mean distortions (solid line) are common while 
maximum distortions to the line profile exceed 20% in most cases. The prob- 
ability to observe a given EW variation was calculated for a similar flow but 
with Cfl ow = 0.7 (shown in dashed line). Variations of the EW by more than 
10% are frequent. 



guish between extrinsic and intrinsic variations of the absorp- 
tion troughs (in the latter case, changes would appear in all 
images after correcting for time-lag effects). 

The results from a systematic study designed to look for 
absorption line variations in lensed quasars may have far- 
reaching implications for understanding the micro-physics of 
the outflowing gas. This would allow us to derive better con- 
straints on the mass loss rate from quasars and, in turn, to 
understand their interaction with the environment. 
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APPENDIX 

Assume that the location and velocity dependent optical depth may be written in the form 

e-^)=T,(R)£(v). (1) 

In this case, 



_ J dvf rfRV (R)/(R) [1 - 77(R)C(v)] _ /■ L g(v)/JR 2 M(R)/(R)T?(R 
JdR 2 fi(R)I(R) J JdR 2 n(R)I(R) 



JdR^(R)I(R) 



(2) 



Clearly, if 77 (R) =const. then the EW does not depend on the magnification fi and is therefore conserved. In case rj is a non- 
trivial function of R then the EW may change as fi evolves. The magnitude of these changes depends on the specific form of 
t(R). For instance, when t(R) is a highly and rapidly variable function of R then effectively non-occulted regions in the disk 
(where t(R) ~ 0) are distributed throughout the disk and are likely to experience different magnifications so that the mean effect 
is suppressed. 



